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A Hamiltonian formulation is given for the gravitational dynamics of two spinning compact bodies 

.4 „„ J ^2 / 4x 



^ , to next-to-leading order (G/c and G /c ) in the spin-orbit interaction. We use a novel approach 

^— ■>) ■ (valid to linear order in the spins), which starts from the second-post- Newtonian metric (in ADM 

' coordinates) generated by two spinless bodies, and computes the next-to-leading order precession, in 

this metric, of suitably redefined "constant-magnitude" 3-dimensional spin vectors Si, S2. We prove 
the Poincare invariance of our Hamiltonian by explicitly constructing ten phase-space generators 
realizing the Poincare algebra. A remarkable feature of our approach is that it allows one to derive 
the orbital equations of motion of spinning binaries to next-to-leading order in spin-orbit coupling 
■ without having to solve Einstein's field equations with a spin-dependent stress tensor. We show 

that our Hamiltonian (orbital and spin) dynamics is equivalent to the dynamics recently obtained 
by Faye, Blanchet, and Buonanno, by solving Einstein's equations in harmonic coordinates. 



PACS numbers: 04.25.-g, 04.25. Nx 



I. INTRODUCTION 



(N 

>. . 

00 , In view of the needs of upcoming gravitational-wave observations, it is crucial to be able to describe in detail the 

■ dynamics of spinning compact binaries. We think that this aim will be fuUfilled by combining the knowledge acquired 

■ by analytical techniques with that obtained by numerical ones. The present paper is devoted to a new, Hamiltonian 
i . analytical treatment of the general relativistic dynamics of spinning binaries. 

T-H • The dynamics of spinning bodies in general relativity is a rather complicated problem which has been the subject of 
many works over many years (starting from the pioneering contributions of Mathisson [J , Papapetrou , Pirani , 
Tulczyjew 14|, and others). This paper focusses on (gravitational) spin-orbit effects, i.e. dynamical effects which are 
linear in the spins of a binary system. The spin-orbit interaction can be analytically obtained as a post-Newtonian 
(PN) expansion. The leading-order contribution of this expansion is proportional to G/c^, while the next- to- leading 
order one contains two sorts of terms: G/c^ and G'^/c^ (here G denotes Newton's gravitational constant and c the 
^ speed of hght). The first complete derivation of leading-order (LO) spin-orbit effects in comparable-mass binary 
' systems is due to Barker and O'Connell [1, @|- These authors derived the spin-orbit interaction by considering the 
quantum scattering amplitude of two spin-^ particles. This curious fact prompted several authors to give purely 
classical derivations of LO spin-orbit effects (see, e.g., Refs. 0, H, [^). For a discussion of LO spin-orbit effects in 
coalescing binary systems see Refs. [lol 

The next-to-leading order (NLO) spin-orbit interaction was analytically tackled only over the last few years. After a 
first incomplete attack due to Tagoshi, Ohashi, and Owen [12], complete results were obtained very recently by Faye, 
Blanchet, and Buonanno [l^ and Blanchet, Buonanno, and Faye [ij. Reference il3] calculated the translational 
equations of motion, as well as the rotational equations of motion for compact spinning binaries to NLO (as here, 
only terms linear in spin were considered). For their derivation, Blanchet et al., working in harmonic coordinates, 
introduced the pole-dipole energy-momentum tensor due to Tulczyjew [3| in the Einstein field equations. They also 
used the general-relativistic-covariant spin supplementary condition (SSC) of Tulczyjew Q or, equivalently in the 
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linear-in-spin approximation, of Pirani [s*] . 

The new derivation of NLO spin-orbit interactions in the present paper is based on a novel approach, and is totally 
independent from the results of Refs. [H, At the end, we shall be able to connect our results to those of [l3l.[l^. 
thereby giving us confidence in the correctness of both investigations. We do not use Tulczyjew's pole-dipole energy- 
momentum tensor. Wc do not cither make use of the Papapetrou (or, more completely, Mathisson-Papapetrou-Pirani) 
translational equations of motion. Our starting point consists of the second post-Newtonian (2PN) metric generated 
by spinless point masses in ADM coordinates, say (7(2pn)o- The crux of our approach then consists in noting that (to 
linear order in the spins) it is enough to compute the NLO spin precession equations in 5(2pn)o to derive the spin- 
orbit NLO contribution in the Hamiltonian, say 7J^'"°(xi, X2, pi, P2, Si, S2). Then, from _ff^'^'-'(xi, X2, pi, P2, Si, S2) 
we can derive the NLO spin-dependent terms in the translational equations of motion (simply by using Hamilton's 
canonical evolution equations). Technically, we shall derive the spin precession equations by starting from the 4- 
dimensional parallel transport equation for the spin 4- vector (with covariant spin supplementary condition) , and then 
by rewriting them in terms of a suitably defined 3-dimensional spin vector, having a constant Euclidean magnitude. 
(This method is essentially that used in Ref. (3] at the LO.) We shall then check the Poincare invariance of our 
Hamiltonian by explicitly constructing ten phase-space generators realizing the Poincare algebra (similarly to the 
proof of the Poincare invariance of the 3PN orbital Hamiltonian given in .15j). After our construction, we shall give 
the relation with the results obtained in Refs. [H, [l^ in the form of explicit transformation formulae. 

We leave to a sequent paper a discussion of the physical consequences of our Hamiltonian formulation, and notably 
its use for improving the description of spin effects within the effective one-body approach [16] . 



II. 3-DIMENSIONAL EUCLIDEAN SPIN VECTOR IN CURVED SPACETIME, AND ITS ANGULAR 

VELOCITY 

When workingto linear order in the spin, the translational and rotational equations of motion of a spinning particle 
in curved space 0, H, 0] (see also |T3| and [11]) read^ 

Du 1 e^^^P ~ 

m—T' ^ o^=SaUpu^R''^;,p, (2.1) 



dr 2 

^ = 0, (2.2) 

where is the normalized 4- velocity of the spinning particle, m^'m^ = —1, m its conserved mass, and S*^ its 4- 
dimensional spin vector; in addition, r denotes the proper time parameter, dx'^/dr = cu^, D the 4-dimensional 
covariant derivative, R^^xp the Riemann curvature tensor, and g the determinant of the 4-dimensional metric 5^,^. 

An important feature of our approach is that we shall not need to consider the translational equations of motion 
(I2.ip . It will be enough to consider the rotational ones (|2.2p . One immediate consequence of (|2.2p is that the 
4-dimensional length of S*^ is preserved along the world line 

g^'^S^S^, = s^ = const, (2.3) 

where g^^ guX — The constant scalar s measures the proper magnitude of the spin. The Eqs. (|2.ip and (|2.2p . to 
linear order in spin, are compatible with the covariant SSC 

S^u" = 0. (2.4) 

At the same approximation, this (Pirani !3]) SSC is equivalent to the Tulczyjew [2] one S^'^p^™ — 0, where = 
mcu^ -\- 0{s'^) is the kinematical momentum (which differs from the canonical momentum we shall use below), and 
where S'^'^ is the antisymmetric spin tensor (see, e.g., [13 )• 

More explicitly, Eq. (|2.2p reads, when expressed in terms of the coordinate time t = a;°/c, 



ds; 

dt 



^ = cT^p'S^v", (2.5) 



^ In this paper, Greek indices run over the numbers 0,1,2,3, Latin indices over 1,2,3, e'^P°'" is the completely antisymmetric (flat 
spacetime) Levi-Civita symbol with e0i23 _ Note that Ref. uses an opposite sign convention for e^f"'^, which leads to an 
opposite sign on the right-hand-side of 1 12.11 1. 
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where F'^^p are the ChristofFel symbols and = c~^dx'^/dt = u^/u^ = (1, w*). Note that, in this paper, we normahze 
the "velocity" = c~^dx^ /dt so that it is dimensionless. 

In addition, we can use Eq. (|2.4p to compute the covariant time component of the spin vector in terms of its 
(covariant) spatial components: 

^0 = ~S,v\ (2.6) 

Substituting this result into Eq. (12. 3|) one finds that the constancy of the 4-dimensional spin magnitude takes the 
3-dimensional form 

G'^S^Sj = (2.7) 

where G"^ is the symmetric matrix: 

G'^ = cf^ ~ g°'v^ ~ g°^v' + g°%'v' . (2.8) 

Now a technically very useful fact is that a positive-definite symmetric matrix such as the one just defined, G*-' , admits 
a unique positive- definite symmetric square root, say H^^ — H^^, such that 

G'^^W^H^^. (2.9) 

This uniqueness result (in some given coordinate system) then naturally leads us to defining a constant-in-magnitude 
3-dimensional Euclidean spin vector Si = 5' as^ 

S, = W^Sj, S,S, = s\ (2.10) 

Upon further use of the spin supplementary condition (|2.6p . the spatial covariant component of the rotational 
equation of motion (|2.5p yields 



dS, 
'dt 

where 



- V^'S,, (2.11) 



;(rJ'.o + rhkv" - r\ov' - r'^.kv'v'') . (2.12) 

Making use of Eqs. (|2.10p and (|2.1ip one can now easily derive an evolution equation for the constant-magnitude 
3-dimensional spin vector Si (dot means differentiation with respect to the coordinate time t): 

Si = V'^Sj, V'^ EE W^{H-^f3 + W^V^\H-^)^^. (2.13) 

The constancy of the Euclidean magnitude of Si implies that the matrix V^^ determining the "rotational velocity" of 
Si = 5' is antisymmetric: V^^ = —V^'^ (a result which is easily checked to hold for the explicit expression of V"^ given 
above). It is then convenient to "dualize" V^^ and to replace it by the 3-dimensional Euclidean (pseudo-) vector 



a = --e^jkV^K (2.14) 



1 

With this notation the rotational equation of motion (|2.13p reads 

S, = + (2.15) 

In other words, we get a Newtonian looking spin precession equation S = f2 x S. 

In summary, the angular velocity of rotation $7 of the constant-magnitude spin 3- vector (|2.10p is directly computable 
from the spacetime metric (and its Christoffcl symbols) by using the explicit formulas (|2.12p . (|2.13p . and (|2.14p . (For 
the self-gravitating spinning particles we are considering, one will need, as usual, to regularize the self-interaction 
terms hidden in the formal results written above. See below.) Note that Jl depends, in general, both on the positions 
and the velocities of all the particles in the system. Indeed, from the explicit formulas above, one sees that Vt depends 
on the velocity of the considered spinning particle. Moreover, the metric and Christoffcl symbols at the location of 
some particle will depend on the positions and velocities of the other particles. 



^ A slightly more geometrical way of phrasing this definition would consist in saying that, starting from a given coordinate system, we 
are constructing a well-defined orthonormal "repere mobile" (or "vierbein") along the worldline of a spinning particle, with respect to 
which the covariant spin 4-vector has components (0, S"). By definition, the spatial components of the metric in this local orthonormal 
frame take the standard Euclidean values &ij , so that we can trivially raise or lower indices on our spin 3- vector. 
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III. DERIVING THE SPIN-ORBIT INTERACTION HAMILTONIAN FROM THE ANGULAR 
VELOCITY OF THE EUCLIDEAN SPIN 3- VECTOR 

Let us now show how the knowledge of the just discussed spin angular velocity vector ft allows one to derive the 
spin-orbit interaction Hamiltonian Hgo, i.e. the part of the Hamiltonian which is linear in the spin variables. 

Let us first recall that a basic result in Hamiltonian dynamics is Darboux's theorem which says that any (non 
singular) symplectic form lu on an even-dimensional manifold can always be (locally) rewritten (after a suitable change 
of phase-space coordinates) in the canonical form uj = "^9^ ^ dp a- When considering N (interacting) spinning 
particles, the dimension of phase space is N{3 + 3 + 2) — 8N, because the description of each particle requires: 3 
spatial coordinates, 3 momenta and 2 spin degrees of freedom, such as two angles 9, 4> needed to parametrize the 
direction of the (constant-magnitude) spin 3- vector S.^. Darboux's theorem then means, in this case, that it is always 
possible to redefine phase-space coordinates such that the symplectic form takes the form 

w = ^ ^ dg* A dpa^ + Sad(- cos Oa) A d(^a) • 

a i 

Here a = 1, . . . ,N labels the various particles (with iV = 2 in our case), while i = 1, 2, 3 labels the spatial dimensions. 
We have written cu in the form it is known to take in special relativity [H, [l^. In the latter case (and, say for 
simplicity, in the case of free particles) , the spin-dependent term in oj was shown to take (globally) the form indicated, 
with Sa denoting the magnitude of the conserved spin of the ath particle, in the sense of (|2.3p . and with 6a and (j)a 
denoting the polar angles of the flat-space limit of the above-introduced constant-magnitude Euclidean spin vector 
Sa, (|2.10p . When considering the interacting case (i.e. turning on a non-zero value of G /c^), and when keeping, for 
simplicity, only the terms linear in spin (so that one can expand the dynamics in powers of both G and Sa), it is 
easily checked (by a perturbation analysis"^) that it is always possible to construct Darboux-type canonical coordinates 
where the spin degrees of freedom are simply the polar angles (in a local orthonormal frame) of the above-introduced 
constant-magnitude Euclidean spin vector S"* 

Finally, we can transcribe this result in the language of Poisson brackets (instead of that of a symplectic form) , by 
stating that there exist phase-space variables x = p = (p°), and S — {Sf) (with a = 1, . . . , A^, and i — 1,2, 3), 

where are, say, the constant-magnitude vectors (|2.10p such that the usual (Newtonian-like) Poisson brackets 

{xl,i^} ^ S^S], {S^5^} = <5"%,5fc^ zero otherwise, (3.1) 

apply to the case of a general-relativistically interacting sytem of N spinning particles. 

Note, however, that this result is essentially kinematical, and has nearly no dynamical content. To describe the 
dynamics of interacting spinning particles, we need to know the expression of the Hamiltonian in terms of the canonical 
variables: H = H{xa, Pa, Sq). As we work linearly in the spins, we look for an Hamiltonian of the general form: 

H(x„ ) i/o(xQ,Pa) +-ffso(xQ, Pa 1 ^a 

). (3.2) 

Here, Ho denotes the orbital part of H, while iJgo contains all the linear-in-spin terms, and can be called the "spin- 
orbit part" . The orbital Hamiltonian Ho is explicitly known up to the 3PN order (isl . [2g | . Our aim here is to compute 
the spin-orbit Hamiltonian Hgo to NLO. Because Hgo is, by definition, linear in the spins we can always write it in 
the general form 

i?so(Xa, Pq, Sa) = ^ i^ai^b, Pb) ' Sq, (3.3) 

a 

where fla — {^D depends on (all) the orbital degrees of freedom (xb,pb), but does not depend on the spins Si,. The 
scalar product in the Eq. p.3p is the usual Euclidean one. 

In Eq. (j3.3p ria is a priori just a notation for the coefficient of in Hgo- But let us now show that it is equal to 
the quantity computed in the previous section, i.e. the angular velocity with which the ath spin vector Sa precesses. 
Indeed, the general principles of Hamiltonian dynamics, together with the canonical Poisson brackets p.ip and the 
form (|3.3p . yield 

Sa = {Sa, Hso{Xb, Pb, St)} = ^aiXb, Pb) X Sa. (3.4) 



^ I.e. by considering general coordinate changes of the form q' = q + 0{s), p' = p + 0{s) and working Unearly in the spins s. 
* As we shall discuss below, we can still modify by a rather general local rotation, but the important point is that our definition of 
S^, I I2.IOI 1. is a smooth deformation of the correct flat-spacetime limit. 
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The only difference between p.4p and the previous result (|2.15|) is that, in p.4p . Jla is expressed in terms of canonical 
positions and momenta, while in (|2.15p it was computed in terms of (say ADM) coordinates and coordinate velocities. 
Because we are working only to linear order in the spin, and because (as was explained above) the canonical phase- 
space coordinates appearing in (|3.3p and (|3.4p differ from the usual ADM-type coordinates used to express the metric 
[and thereby to compute the angular velocity r2a(x^^^, v^^'^) by means of (|2.12p . (|2.13p . and (|2.15p ] only by terms 
proportional to the spins, it suffices to use the known [l^llol spinless link between ADM momenta and ADM velocities 
to compute na{xb,Pb) from flaix^^^ ,v^^^). 

In the previous section we introduced a specific, well-defined "conserved" spin 3-vector Si to parametrize the two 
degrees of freedom of a spinning particle. Our choice had the nice features of being universally associated to the 
choice of a coordinate system, and of reducing to the choice made in the flat spacetime limit fT9'|. However, it was by 
no means physically unique. 

Let us now show that the freedom in the choice of conserved spin vector is simply a "gauge freedom" (local rotation 
group) which does not change the physical results one can deduce from the Hamiltonian. Indeed, the condition 
SiSi = leaves as ambiguity in the definition of the conserved spin variable Si a local 3-dimensional Euclidean 
rotation Si S^, with 

Sj^ — RijSj^ (*^-^) 
where R is an arbitrary rotation matrix. It is sufficient to consider the case of an infinitesimal rotation, say 

Rij — Sij — 9ij, (3.6) 
where 6ij is a small antisymmetric matrix. This leads to an infinitesimal change 

6S = exS, (3.7) 

where we introduced the dual vector 6 such that 6ij — SijkOk- 

Let us show that such a change can be considered as being induced by an infinitesimal canonical transformation g 
in the full phase space (x, p, S). (Canonical transformations are symmetries of Hamiltonian dynamics. In particular 
they preserve the basic Poisson brackets written above.) We recall that such a canonical transformation acts on any 
phase-space function / according to 

5f^{f,9}. (3.8) 
It is then easily checked that a transformation of the form 

5(x,p,S) = 0(x,p) -S (3.9) 

transforms the spin vector according to 

(5S = {S,g} = 6/ X S, (3.10) 

which exactly reproduces the effect of an infinitesimal local rotation written above. However, we have learned that 
such a local rotation must be accompanied by a corresponding transformation of the orbital degrees of freedom (x, p) 
of the form: (5x — {x,(7}, i5p = {p,^}. Then, under the simultaneous changes of x, p, S induced by the canonical 
transformation g (and the corresponding change of the spin angular velocity 17' ~ + A9/dt) one finds that the 
numerical value (evaluated at corresponding phase-space points) of the Hamiltonian is invariant. 

We have therefore shown that the arbitrariness in the "rotational state" of the conserved spin is simply (as expected) 
a "gauge symmetry" (under a local S0(3) group). 



IV. DERIVATION OF THE SPIN-ORBIT HAMILTONIAN IN ADM COORDINATES 



Let us now sketch the computation of the NLO angular velocity fii in ADM coordinates [which will then give us 
the NLO spin-orbit Hamiltonian according to Eq. p.3p ]. 

As usual we split the four-dimcnsional metric into three-dimensional objects {a, Pi,"fij), where 

a=(-g"0)-i/2, (3, = go,, i,,=g,j. (4.1) 

One can show, using the definitions /3' — 7*-'/3j , 7'-'7j7c = <^ife, that the following exact formulas hold: 



-^0 



i Oi 

a 



+ (4.2a) 
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r°„ = ixy, (4.2b) 

= ^ySfe^.o - l^P'aj + \l'\Pk,j - f3j^k) - l^fi^P^Kkj, (4.2c) 

r^k^'r),--(3^K,k, (4.2d) 

a 

where Kij is the extrinsic curvature of the constant time shce. Note that, for convenience, we use the Kij ^ + jlj 
sign convention (instead of the — 7ij convention used e.g. in Ref. [l7j). In terms of the field momenta tt'^ it reads, 



IGttG 1 / 



V7 



{j^kl3l-llvlkl)7T''\ (4.3) 



where 7 = det(7y). The Christoffel symbols related with the 3- metric 7ij are denoted by r*j.. Let us also note 
that the dimensionless coordinate velocity can be expressed in terms of the bare kinematical linear momenta 
pbarc _ j^gy.^ generality, as follows 



ij ^baro 

Note, however, that the latter result applies to the canonical momentum only modulo corrections proportional to the 
spin. 

We employ the ADMTT (ADM transverse-traceless) coordinate conditions [2l| 
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7^, = (^l+g'^j S,,+ht/, 7r"=0, (4.5) 
and recall that 

with TT^rp being of the order 1 /c^ . 

Let us now expand all quantities in a post-Newtonian (PN) expansion. Here and below the subscript (n) indicates 
the part of a quantity which is of the nth post-Newtonian order, i.e. which is proportional to (1/c^)". For instance 
we decompose: 

Here ri(2)i c>c G/c^ is the well-known LO contribution 0,0, BIl], while ri(4)i oc G/c'^ -f-G^/c^ is the NLO contribution 
that we wish to compute. These contributions arc more explicitly given in terms of the "precession velocity" V^-' of 
the "coordinate spin vector" Si, which entered Eq. (|2.1ip . Inserting in Eq. (|2.12p the Christoffel symbols (|4.2I) . and 
then inserting the result in Eq. (|TTn)l [where H'-i is computed from Eqs. ([^ . ([^ . (gT]), (g^l), (^3)) . and (g^], we 
obtain the following more explicit formulas for the 3- vectors fl(^2)i and 51(4)^ from Eq. (|4.7p : 

^^(2)z/c = ^%fc (^/3(3)j,fe + (a(2)j - ^0(2) , (4.8a) 

„.k lTT 

^(4)i/c - -%fc ( /5(5)j,fe + /?(3)fca(2)j - 2<^(2)/3(3)j,fe + Y^0(2)'^(2)j^' " -j^<t>(i),j^ " ^(4)fci,j" 

+ (a(4) J - a(2)a(2),j)i;'' -f ^^g)"''"' - ^a(2),fct'^w'w' + J^w'^w'i;'^ . (4.8b) 

At this point, it only remains to implement three technical steps: (i) to insert the explicit form of the 2PN-accurate 
metric describing two spin- less particles in ADMTT coordinates (from [2^ and [13]), (ii) to replace the velocities by 
their IPN-accurate expression in terms of the canonical momenta pi, and, finally, (iii) to regularize the self-interaction 
terms that arise when evaluating Eqs. (|4.8p . 

The explicit expressions for the metric functions entering Eqs. (|4.8p can be found e.g. in Appendix A of Ref. [2^ 
(where the functions 0(2), '/'(4), '''(3)7 ^.nd hj^ij can be found) and in Ref. 23] (where the functions Q!(2), Q!(4) and 
/3(3)i, /3(5)i are given). 
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As for reexpressing the velocities in terms of raomenta, it yields a further PN-expansion of the form = ^'^(i) + 
^a(3) where^ '^a{i) coordinate velocity of the ath particle expressed in terms of the canonical variables 

Xa and Pa at the Newtonian accuracy, i.e., i^^f^i-j — Pai/ {m.ac), and w^j-g-j is the IPN correction to Pai/ (mac). The 
latter IPN correction explicitly reads 

CfniT-po) ; / P? SGmo \ 7G 



the expression for V2{'i) can be obtained from the above by exchanging the particles' labels. 

The final step then consists in evaluating (note that the meaning of ^^^^2) ^nd ^^^(4) slightly different because 

of the re-expansion of velocities in a PN expansion) 

K{2)/C = ^Ey-fcRega (^/3(3)j- fc + (a(2)j - ^0(2)j)t'a(l)). (4.10a) 
^a(4)/c= ^£yfcR.ega(^A5)i,fc +/3(3)fca(2)j " ^'/'(2) /3(3)i,fe + ^0(2)0(2) jWa(l) " ^'/'(4),i^'a(l) " ^(4Tfci,i^^i(l) 



+ (a(4)j - a(2)a(2)j)Wa(l) +*(3)^a(l)^'i(l) - 2"(2):fe"a(l)«i(l)^'i(l) + 4 ^"a(l)«i(l)«i(l) 

+ («(2),.-^'^(2),,)<3)), (4-lOb) 

where Rega(/(x)) indicates that one must regularize the limit x — > x^. At the level at which we are working, 
this regularization is not ambiguous and can, for instance, be simply performed by using Hadamard's "partie finie" 
regularization (as explained, e.g., in Appendix B of Ref. ^^)- The final results we got read 

_ G ( 3m2 \ 

iii(2) = -^-^ 7^ ni2 X Pi - 2ni2 X p2 I , (4.11a) 

G2 / / 11 ( 15 



^^1(4) = ^I^r - y ™2 - 5— j ni2 X pi + (^6mi + —m2 j ni2 x p2 

G If 5m2pj 3(pi-p2) 3p^ 3(ni2 • pi)(ni2 • P2) 3(ni2 • P2)" 



„ s 9 ' 9 M^12 ^ Pi 

'^^12 V V Sm'i Amf 4mim2 4to^ 2TOirn2 / 

(Pi-P2) , 3(ni2 • Pi)(ni2 • P2)\ , /3(ni2-pi) 2(ni2-p2)\ \ 

\ ni2 X p2 + Pi X p2 . (4.11b) 

mim2 TOir7i2 / \ 4rnj TO1TO2 / / 

The expressions for r22(2) and r22(4) can be obtained from the above formulas by exchanging the particles' labels. 

From these results we can then explicitly write the spin-orbit Hamiltonian to leading and next-to-leading PN orders. 
Indeed, 

Hso{^a,Pa,Sa) = ^^a{xb,Pb) ' = ^ {^a{2){^b , Pb) + ^a{i){^b,Pb)) ' S^. (4.12) 
a a 

More explicitly, the separate LO and NLO contributions in the PN expansion of the spin-orbit interaction term, 



iJ,o(Xa,Pa,Sa) = li/LO(x^^p^,Sa) + ^ H^J''' i^a,Pa,Sa) + O ] , (4.13) 



^ Here and below a,b = 1,2 are the particles' labels, so ma, Xa = (x^), and pa = (pai) denote, respectively, the mass parameter, the 
position vector, and the linear momentum vector of the ath body; for a ^ b we also define r^^ = Xa — X;,, r^t = \j^ab\, ^ab = ^ab/^abi 
I ■ I stands here for the Euclidean length of a 3-vector. 
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read, 

Hlfi^a,Pa, Sa)=c'Y, ^a(2)i^b, Pb) ' S,, (4.14a) 

a 

i/f,L°(x,, p., Sa)=C^Yl ^alM^b, Pb) ■ Sa. (4.14b) 

a 

Finally, note a remarkable feature of our Hamiltonian approach to spin-orbit effects: the sole computation of the 
rotational velocity of the ("conserved") spin vector (given by parallel transport in the 2PN-accurate metric of N 
spinless bodies) determines the NLO spin- dependent terms in the translational equations of motion of N spinning 
particles. Indeed, the sole knowledge of Qa{^b,Pb) yields that of the total spin-dependent Hamiltonian p.2p with 
p.3|) , so that the general principles of Hamiltonian dynamics (with canonical Poisson brackets) yield 

dH{xb,Pb,Sb) . ai/(xb,pfc,Sb) 
= + d^. ' = d^. • ^'-''^ 

In view of the availability of algebraic manipulation programmes, there is no need to write down explicitly the trans- 
lational equations of motion (|4.15p . with NLO accuracy in spin-orbit terms (and 3PN accuracy in spin- independent 
terms [H,!!^). We shall verify below that the Hamiltonian, ADM-coordinate translational equations of motion (|4.15p 
are equivalent to the harmonic-coordinate ones recently derived in (isl . [l3 | by a more complex calculation which 
involved the explicit consideration of spin-dependent contributions in the metric. 

V. POINCARE INVARIANCE 

The general relativistic dynamics of an isolated A^-body system should admit the full Poincarc group as a global 
symmetry (because it is a symmetry which preserves asymptotic flatness). On the other hand, this symmetry is 
not manifest in the Hamiltonian ADM approach to the A'^-body dynamics because it splits space and time, and uses 
non-Lorentz covariant coordinate conditions. In a previous paper [l5|, treating non-spinning particles, the authors 
showed how to bypass this technical mismatch: the basic idea is that, in the Hamiltonian formalism, the global 
Poincare symmetry is realized in phase-space in a non-linear manner. However, one can efficiently detect the presence 
of this symmetry by proving the existence of ten phase-space generators i?(xa, Pa, Sa), Pi(xa, p^, Sa), Ji(xa, Po, Sa), 
Gi(xa,Pa,Sa) (depending on all phase-space variables) whose Poisson brackets reproduce the standard Poincare 
algebra. In the case of non-spinning particles, Ref. constructed the ten generators of the Poincare group at the 
3PN level of approximation. We shall show here how to extend this construction to the more involved case of a system 
of spinning particles. 

Let us first recall the explicit Poisson-bracket form of the Poincare algebra that should be realized: 

{P„Pj} = 0, {P,,H}^0, {J„if} = 0, (5.1a) 
{Ji,Pj} = EijkPk, {Ji,Jj}=SijkJk, (5.1b) 
{Ji,Gj} ^ EijkGk, (5.1c) 
{G„i7}-P„ (5.1d) 

{G,,Pj} = ^HS,,, (5.1e) 

{G„Gj} = -^£yfeJfe. (5.1f) 

The translation. Pi, and rotation, J^, generators are simply realized as 

Pj(Xa,Pa,Sa) = ^Pai, (5.2a) 

a 

Jj(Xa,Pa,Sa) = ^ (siklx'^Pai + Sai) ■ (5.2b) 
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Note the very simple, additive, form of these generators, and, in particular, how our Hamiltonian "conserved spin" 
variables appear as Newtonian- like (but relativistically correct) contributions. 

As for the Hamiltonian we already know that (in our linear- in-spin approximation), it is a sum of an orbital 
part. Ho, and of the above-determined spin-orbit part, Hgo, Eqs. (|4.12p . (|4.13p . and (|4.14p : 

Pa; 

). (5.3) 

The orbital Hamiltonian Ho (including the rest-mass contribution) is explicitly known up to the 3PN order p^ . [20j : 

Ho{Xa,Pa) = ^maC^ + HoN{Xa,Pa) + ^ -ffolPN (Xa , Pa) + ^ -ffo2PN (Xq , Pa) + ^ i?o3PN (Xa, Pa) +0 ( . (5.4) 

The most delicate generator to consider is the boost (or center-of-mass) vector G. It can be represented as a sum 
of "orbital" and "spin-orbit" parts 

G(Xa, Pq, Sa) = Go(Xa, Pa) + Gso(Xa, Pa, Sa), (5.5) 

where, as everywhere in this paper, we call "spin-orbit" the part which is linear in the spin variables. The orbital 
part, Go, was explicitly determined up to the the 3PN order in Ref. [l5|: 

Go(Xa, Pa) = ^ m-aXa + \ GolPN(Xa, Pa) + \ Go2PN(Xa, Pa) + \ Go3PN(Xa, Pa) + O [ ^ j . (5.6) 



The spin-orbit part can be decomposed in leading-order (LO), next-to-leading-order (NLO), and further contributions: 



G,o(Xa, Pa, Sa) = ^ {^a, Pa, Sa) + ^ G^o"" {^a, Pa, Sa) +o[^\. (5.7) 



The leading-order term in (15.71) is known from the special-relativistic limit (by replacing the special-relativistic energy 
in the results of, e.g., Refs. [l9l[25j. by the rest-mass contribution) 

G^O(xa,Pa,Sa) = - + (1 ^ 2), (5.8) 

Lm\ 

where the operation "-|-(1 ^ 2)" denotes the addition to each displayed term of another one obtained by exchanging 
the particles' labels. 

The real difficulty lies in constructing the NLO contribution to the boost generator (and in proving that it satisfies 
the correct Poincare algebra displayed above). We solved this problem by using (as in our previous work [15]) 
the method of undetermined coefficients. The most general form of G^'^'-' can a priori depend on eight unknown 
dimensionless numerical coefficients 171, . . . , 

pNLO _ P? a „ „ , / Si X pi Si X pz / (ni2 -Pi) (ni2 ■P2) \ 

*^so — T, — 3^1 X Pi H 51 r 52 V gz V 34 ni2 X hi 

oTOf ri2 \ mi 7712 \ mi m2 / 

/ (5i,ni2,_pi) (51,7112, P2)\ \ , Gm2 / (5*1,7112,^1) (^i, 7112,^2) \ , 
' ' 56^ ^ ni2 -F^— 57^ -"Vg%- ^ xi + (1^2), (5.9) 



\ TOi m2 / / r\2 \ mi TO2 

where we have introduced the following notation for the Euclidean mixed product of 3-vectors: (Vi , V2 , V3) = Vi • (V2 x 
V3) = ^ijkViV^V^ . Note that the coefficient of the first term on the right-hand-side is determined by considering 
the special-relativistic limit [l^, [l^, [2^ . We have also used some structural information coming from a conceivable 
field-theory computation of G (say as the space integral of the Oi component of some effective stress-energy tensor). 
Indeed, such a computation could be thought of in terms of some Feynman-like diagrams, where the interaction terms 
(i.e. those containing a power of G) would all be proportional to some basic "source" term involving either Si and 
m2 (connected by a propagator, and possibly some power of the velocities, vi ^ pi/mi or V2 ~ P2/m2), or similar 
terms involving 52 and toi. The main point being that pure "self-interaction" terms (say proportional to Si and toi) 
cannot appear. 

Let us now consider the explicit Poincare algebra requirements of Eqs. (|5.1ap - (|5.1f|) . It is easily verified that the 
generators Pi, Ji, H, and Gi, in the forms given above, exactly satisfy the relations (|5.1ap . (|5.1bp . and (|5.1cp . We now 
consider whether the center-of-mass vector G with the 2PN spin-orbit part given by Eq. (|5.9p can satisfy the three 
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relations (|5.1dp - (j5.1f|l . This requirement yields many equations that have to be satisfied by the unknown coefficients 
gi, . . . ,gg,. We have first found that there exist unique values of the coefficients gi, . . . ,gs ensuring the fulfillment of 
the sole relation (jS.ldp . These values are 

5 3 11 3 

91 = |, 52 = 93 = 0, 34 = --, 55 56 = 1, 97 ^ 98 = -2- (5.10) 

Then we have checked that the solution (|5.10p also guarantees the fulfillment of the remaining relations (|5.1ep and 

dm. 

In summary, we succeeded in proving the Poincare invariancc of the above-defined NLO spin-orbit interaction 
[determined by Eqs. (I4.12p . (|4.13p . and (|4.14p ] by explicitly constructing ten phase-space generators satisfying the 
Poincare algebra brackets of Eqs. (|5.1ap - (|5.1f|l . 



VI. COMPARISON WITH HARMONIC-COORDINATE-BASED RESULTS 

References [l3l [l^ recently computed, by means of two separate calculations and in harmonic coordinates, both 
the NLO spin-dependent contributions in the translational equations of motion of two spinning particles, and the 
corresponding NLO terms in the spin precessional equations of motion. In the present Section we shall prove that 
our results are physically equivalent to the results of Refs. [l^, [13] by finding the explicit form of the transformation 
that match the ADM variables used by us with the harmonic variables used in Refs. (T3l.[T^. Let us start by warning 



the reader that in the whole paper [13| and in most of the paper [lj| Blanchet et al. chose to express their results 
in terms of some "non-conserved" spin variables S]^^^, i.e. variables whose Euclidean magnitudes are not conserved 
in time. It is only in Sec. VII of that redefined spin variables with conserved Euclidean lengths, say SJj^^^, are 
introduced and used. 

Our task here will be to exhibit the explicit transformation between the "ADM variables" (xa,Pa,Sa) used in our 
work, and the "harmonic variables" (ya, Va = ya, Sl^^^) used in [ll,[l3], and to prove that this transformation maps 
the two sets of results into each other. (It is more convenient for us to exhibit the link with the "conserved" version 
of the harmonic spin variable used by Blanchet et al. The relation between their two spin variables, S^^^ and SJij^^^, 
is given in Eq. (7.4) of [l4|.) 

We write the transformation of variables in the general form^ 

yait) - YaiMt), Pbit), Sbit)), (6.1a) 

ScBBF(^) ^ S,(xf,(t), pb(t), Sf,(t)). (6.1b) 

Let us first find the transformation between spin variables. Section VII of Ref. [l3| gives (see Eq. (7.6) there) 
the explicit result for the angular velocity vector ^2^^^ of their conserved harmonic spin variable, yielding a spin 
precessional equation of motion of the form 

jQcBBF 

^^^ = n^^^y<S:^^^, a =1,2. (6.2) 
at 

They give the NLO expression of fJ^^^ in terms of the harmonic orbital coordinates (ya, v^). We have re-expressed 
r2j^^^(yh, Vf,) in terms of ADM coordinates and momenta, to IPN accuracy (using the well-known link between the 
two sets of variables^). We then compared the result with our results (|4.1ip . We have found 

^a(2^iyb,^b) = rJa(2)(xb,Ph), (6.3a) 

d6)„ 



^ai4 (yfc, Vb) = f^a(4)(xb,Pfc) + (6.3b) 



where 



G f (ni2-p2) , (ni2-p2) 9 > 
di = —. ni2 X pi H ni2 X p2 - pi x p2 . (6.4) 

c*ri2 \ 4mi m2 4mi ' 



° Here, both sides refer to the same numerical value of their respective coordinate times. 

^ We recall that harmonic and ADM coordinates coincide at IPN, but that one must transform velocities into momenta by means of the 
IPN transformation Eq. 114.91 1. See [24{ for the 3PN-accurate version of this transformation. 
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From the results (|6.3p -- (|6.4p it is easy to deduce that the two sets of spin precession equations of motion are physically 
equivalent* , and that the two sets of spin variables are related as in the general transformation links written above 
with a spin transformation Ha of the explicit form: 

Sa(xb, Pb, Sf,) =Sa + 0Q(xb, Pb) X Sa- (6.5) 

In other words, our conserved spin variable differs from the conserved spin variable defined in Eq. (7.4) of J^] by a 
small (time-dependent) rotation of angle 9a{^b,Pb)- Such a difference was a priori to be expected because constant- 
magnitude spin vectors are not uniquely defined. We have shown above that to each choice of coordinate system is 
canonically associated a particular choice of local orthonormal frame (along the worldline of a spinning particle) , and 
thereby a particular choice of "conserved" spin 3-vector. We have investigated whether the conserved spin 3-vector 
defined by Blanchet et al. does correspond to applying our general definition to the case of harmonic coordinates. The 
answer is "no" . We found that if Blanchet et al. had used our general definition ()2.10p in their harmonic coordinate 
system, the angular velocity fla that they would have obtained would differ from our ADM spin vector by a rotation 
vector 6a differing from the result above by having the factor 9 replaced by 1 in the last term of Eq. (|6.4p . There 
is nothing surprising in such a difference as the spin re-definition used by Blanchet et al. was somewhat arbitrary. 
Anyway, as already mentioned above physical results will not depend on such "gauge choices." 

Let us now turn to the determination of the transformation Yq between ADM and harmonic orbital degrees of 
freedom. As usual we can decompose Yq into spin- independent, Y°, and spin-dependent (and linear-in-spin), Y^°, 
terms: 

Ya(xb,Pb,Sb) -Xa+Y°(xb,Pb)+Y^°(xb,pb,Sb), (6.6) 

where the spin-dependent term is of the form 

Y-(xb, Pb, Sb) - Y:'^^) i^b, Pb, Sb) + Y^-^^) (xb, Pb, Sb) + 0(c-6). (6.7) 

The spin-independent part of the transformation was explicitly given, up to the 3PN order, in Ref. The leading 
order spin-dependent part has been known for many years (see, e.g., Ref. [2^). and equals 

Y:^,)(xb,Pb,Sb) = %^. (6.8) 

We have determined the next-to-leading order spin-dependent part, Y^°^^ , by using again the method of undetermined 
coefficients. We have considered the most general template for Y^°^^ which depends (after using the special relativistic 

limit to determine the 1/c* term which remains in the G limit, and structural information of the same type as 
that explained above in the case of Ggo)^ on 12 unknown coefficients. It reads 

vso / c ^ Pi c , ^"^2 ^ ( Si X pi Si X p2 

Yim(Xa,Pa,Sa) = , 4 S1 Xpi + — ai -1-02 

^ ' OC*mJ C 7^12 TTll \ mi 7712 

(ni2-pi) ^ (ni2-p2)\ „ , ( {Si,ni2,pi) {Si,ni2,P2)\ \ 
as- - + aA- ni2 X Si -I- 05^ + ae- ni2 



nil 1712 J \ mi 1712 

G / S2 X pi S2 X p2 / (ni2-pi) (ni2-p2) 

Ol 1- O2 h O3 h O4 

c*ri2 \ mi m2 V mi TO2 



, {S2,ni2,Pl) (52,ni2,P2) 'i , , , 

+ [b5- ^+^6^ ^ ni2 . (6.9) 



mi m2 



One can now think of two different ways of determining whether there exists a set of coefficients ai, . . . , ag; 61, . . . , 65 
such that our translational Hamiltonian equations of motion (with NLO spin-dependent terms), Eq. (|4.15p . are 



* As a further check, we have also explicitly verified that the Hamiltonian time derivative (computed with our dynamics, namely 
{Sj^^^,-ff}) of the originally defined (non-conserved) spin vector Sj^^'^ of [T3II coincides with the NLO spin precession law given 
by Eqs. (6.1)-(6.3) there. To do this calculation we defined the phase-space quantity Sj^^^(x(,, p;,, S;,) by inserting Eqs. 1 16.31 1. 1 16.41 1 into 
Eq. (7.6) of [ij. 

® More specifically we required that, say, miY|° be proportional (modulo some velocity-dependent factors involving Va ~ Pa/ma) either 
to m\S2 or to m2S\. 
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physically equivalent to the corresponding translational harmonic equations of motion derived in [1J|. (1) A first way 
would consist of inserting the putative general transformation Yo(ai, . . . , ag; 61, . . . , fog) directly into the translational 
equations of motion derived in [l^] (using the fact that we have already determined how their spin variables are linked 
to ours), and to compare the result to the explicit form of our translational Hamiltonian equations of motion, Eq. 
(|4.15[) . This approach is, however, computationally heavy. (2) Therefore, we have instead used a simpler approach 
consisting in comparing the ten conserved quantities derived in harmonic coordinates in Ref. 13], namely the energy 
E{ya,^a,^^^^), the total linear momentum P(ya, Va, S^^^), the total angular momentum J(ya, Vq, Sj^^^), and 
the center-of-mass vector G(ya, v^, S)^^^), with the ten phase-space Poincare generators constructed above within 
our Hamiltonian formalism. To do this comparison explicitly, we first need to perform two replacements: (i) to 
replace the non-conserved spin variable S^^^ used in [l3| in terms of the conserved one SJ^^^^ introduced in |14|, 
thereby obtaining new expressions E{ya,Va, 8^^^^), P(ya, v„, S^^bF)^ j(y^^^ g^.^BF)^ Q(y^^ S^bbF) for the ten 
conserved quantities, and then (ii) to replace the harmonic-coordinate velocities = dya/di in terms of Hamiltonian 
time-derivatives, namely ~ {Ya, H}- Finally, the values of the coefficients oi, . . . , and 61, . . . , must fulfill 
the equations 

E{Ya{Xb, Pb, Sb), Va(x6, Pb, Sb), Sa(xb, Pb, Sfc)) = (x^, p^, Sq), (6.10a) 

P(Y,(xb, Pb, Sb), V„(xb, Pb, Sb), S,(xb, Pb, Sb)) - 51 P-^' C^.lOb) 

a 

j(Ya(xb,Pb,Sb), Va(xb,Pb,Sb),S<j(Xb,Pb,Sb)) = ^ (Xa X p<j + S^) , (6.10c) 

a 

G(Ya(xb, Pb, Sb), Va(xb, Pb, Sb), Sa(xb, Pb, Sb)) = G(Xa, Pa, Sq). (6.10d) 

By considering the first three of these equations (i.e. by comparing the two expressions for the energy, the total linear 
momentum, and the total angular momentum), we obtained a unique set of values for all the unknown coefficients 
fli, . . . , ag; 6i, . . . , feg. We then verified that these values satisfy also the fourth of Eqs. (|6.10p (thereby giving us 
confidence in the correctness of our Hamiltonian, and providing many non-trivial checks of the previous results 

MM)- 

Our unique solution for the spin-dependent transformation of orbital coordinates Y^°2) + ^a°4) rG^-ds: 

Si X pi Si X pi / p^ Gni2\ 



Yi'(2)(xQ,Pa,Sa) 4- Y^4)(Xq,Pq,Sq) - 2c^ml '^^^ 

G 



1 d^m\ \ 8mi ri2 



2c* 



m2ri2 



3S2 X p2 + 2(ni2 • P2)ni2 X S2 -f (5*2, ni2,P2)ni2 ) . (6.11) 



Note that the first three terms on the right side of Eq. (|6.1ip (i.e. the terms proportional to Si) have the same 
structure as the exact special relativistic value [Tsl [Tol. [25j for the shift Yf° between the canonical^° orbital coordinate 
xi and the usual Lorentz-covariant (harmonic) orbital coordinate yi , namely 

yi=xi + „ . , (6.12 

TOi(miC'^ -I- i/i) 



where Ei = -(/(mic^)^ -I- (pic)^ is the relativistic energy (including the rest-mass contribution). The pi-dcpendcnt 
terms in the first three terms of Eq. (I6.1ip correspond to the NLO expansion of the special-relativistic result, while the 
additional G-dependent contribution can be roughly understood as a gravitational addition to the special-relativistic 
energy Ei (though it does not have the correct coefficient to be really interpreted so simply). 

By contrast, the terms proportional to S2 in Eq. (|6.11|) do not have correspondants in the special-relativistic (i.e. 
G — !■ 0) limit. As was to be expected they vanish in the limit where the second body (of mass 7712) is heavy and 
fixed (p2/77i2 — > 0). (Indeed, if we consider a non-spinning test particle, mi, 5*1 = 0, moving in the background of 
a fixed, heavy spinning mass, m2,S2, the harmonic-coordinate geodesic action of mi will already yield a canonical 
Hamiltonian action.) We leave to future work a direct derivation of these terms from the perturbative construction 
of canonical coordinates (of the type q — y + 0{s), p = + 0(s)) alluded to above. 



The classical canonical variables, here denoted Xq, P(,,Sa, correspond, at the quantum level, to the so-called Pryce-Newton-Wigner 
variables. 



13 



Finally, as a further check on the algebra, we have also used the "direct" method (1) mentioned above (the first 
method we could have used to determine the values of the coefficients Oi, . . . , ag; &i, . . . , 65). More explicitly, we started 
from the harmonic-coordinate translational equations of motion with NLO spin-orbit effects given in Eqs. (5.3) of Ref. 
|13|. We then replaced in these equations the non-conserved spin vector S^^^ by its expression (as given in Eq. (7.4) 
of [ij]) in terms of their conserved spin vector S^^^^. This yields 2PN-accurate translational equations of motion of 
the form 

^ = AN^(y,,v,) + l(Air(y.,v,) + A^o (y,,v,,SgBBF)j 

+ ^(Ar„^(y.,v,)+Afr(y.,v,,SgBBF))+o(e-6). (6.13) 

We then compared the right-hand-side of Eq. (|6.13p . let us denote it by Aa — Aa(y6, v;,, S|;^^^), to its direct 
Hamiltonian recomputation by means of our Hamiltonian flow, i.e. 

A, = {V„, H} = {{Y„ H}, H), (6.14) 

together with the needed transformations ()6.ip (determined above) between harmonic and canonical variables. Again, 
this verification worked perfectly and (together with the similar direct verification of the NLO spin precession equation 
mentioned above) gives us confidence that both sets of results (harmonic and Hamiltonian) are correct. 
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